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ABSTRACT 

Evolution of first population of massive metal-free binary stars is followed. Due to the low metallicity, 
the stars are allowed to form with large initial masses and to evolve without significant mass loss. 
Evolution at zero metallicity, therefore, may lead to the formation of massive remnants. In particular, 
black holes of intermediate-mass (~ 100 — 500 M-q) are expected to have formed in early Universe, in 
contrast to the much lower mass stellar black holes (~ 10 M.q) being formed at present. Following 
a natural assumption, that some of these Population III stars have formed in binaries, the physical 
properties of first stellar binary black holes are presented. We find that a significant fraction of such 
binary black holes coalesces within the Hubble time. We point out that burst of gravitational waves 
from the final coalescences and the following ringdown of these binary black hole mergers can be 
observed in the interferometric detectors. We estimate that advanced LIGO detection rate of such 
mergers is at least several events per year with high signal to noise ratio (> 10). 
Subject headings: binaries: close — black hole physics — gravitational waves 



1. INTRODUCTION 

The properties of Population HI stars have stirred a 
lot of interest in the recent years. It has been real- 
ized that zero metallicity stars with mass es up to sev- 
eral hui idred of solar ma. s ses ar e stable ijBaraffe et al.l 
12001'). 'H eger fc WooslevI lj2002D estimated black hole 
(BH) masses formed in the evolution of high- mass metal- 
free stars. For initial masses above 40 Mq the remnant 
is a BH with the mass essentially the same as the pro- 
genitor with the exception of stars with initial masses 
between 140 and 260 M.q which undergo a pair insta- 
bility s upernova (SN) exp losions and leave no remnant 
at all. iHeeer et al.l l)2003() found that these conclusions 
hold also for low, non-zero metallicity stars. Numer- 
ical studies of collapse and fragmentation of a metal 
free gas in the early Universe ijBromm 6^1111999112001 
[Omukai & Palla 2003) indicate that the initial mass func- 
tion (IMF) of Population III stars is top heavy, and might 
be bimodal w i th the high mass peak aroun d 100 
llLarsonI Il998t iNakamura fc Umemural 120011: iChabried 
l200,3j) . Everv known stellar system contains a consider- 
able fraction of binaries. Therefore it seems quite natural 
to allow for a possibility that some fraction of massive 
Population HI stars formed in binaries as well, and then 
investigate their subsequent evolution. There is a possi- 
bility that these systems form massive black hole black 
hole (BH-BH) binaries which in turn may be observable 
in gravitational waves by the interferometric detectors. 

In this paper we consider the properties of Popula- 
tion HI binaries with the initial component masses in 
the range 100 — 500 Mq. In § 2 we describe the model of 
evolution of metal-free systems leading to formation of 
massive BH-BH binaries. In § 3 we present characteristic 
properties of the BH-BH population and discuss the ob- 
servability of their mergers. Conclusions and summary 
of results are given in § 4. 



2. EVOLUTIONARY MODEL 

A simple model of the evolution of Population III stars 
was constructed. Lacking the observational input, we 
used a set of recent calculations for metal-free single stars 
and then combined them with the basic binary evolution- 
ary prescriptions. 

Using the numerica l calculations of ste l lar tracks for 
Popu l ations HI stars jBaraffe et alJl2nnH iMarigo et alJ 
1200 It iSchaerej l2002|) we obtained the evolutionary 
timescale and radial expansion history of a star as a func- 
tion of its initial mass. After the core hydrogen exhaus- 
tion, we calculate He-core core mass fo r a given star usin; 
the ap proximate empirical formula of iHeger fc Woosle 
l)2002t) . Once a sta r has finished its nuclear evolution, we 
follow the work of IHeger et al.l (|2003) to decide on core 
collapse outcome. Depending on the initial star mass, the 
low-metallicity massive star may either collapse directly 
and form a BH (without accompanying SN explosion) or 
be entirely disrupted (no remnant left) in a pair instabil- 
ity SN (stars within initial mass range of 140-260 M©). 
In the former case, the total mass of collapsing star forms 
a BH, and we assume that there is no natal kick asso- 
ciated with the direct BH formation. The pulsational 
pair instability is also taken into account as it may re- 
move part of the s tar envelope just prior to the collapse 
(|Heger et al.l 12001 . We assume that half of the enve- 
lope is lost for the stars with initial mass in the range: 
100-140 Mq. 

The orbit of each binary is assumed to be circular, and 
the orbital separation is drawn from a distribution flat 
in log arithm with the maximum value of 10^ R© l|Abd 
Il983j) . Since little is known about the shape of the IMF 
of Population HI stars we adopt a power law shape with 
a = —2 exponent and draw the primary (initially more 
massive component) mass from such a distribution in the 
range of 100 — 500 M©. The secondary mass is obtained 



2 



Belczynski, Bulik & Rudak 



BH- 



Channcl 

bhbhOl 
bhbh02 
bhbhOS 
bhbh04 
bhbhOS 

^BH1/BH2: first/second BH formation. CE: common enve- 
lope, MT: non-conservative RLOF, where 1 or 2 denotes the 
donor, either primary or secondary, respectively. 



through mass ratio (secondary/primary) which is drawn 
from a flat distribution. Radial expansion of the compo- 
nents is foUowed, and in the case of Roche lobe overflow 
(RLOF) we apply one of the following prescriptions to 
calculate the outcome. For unevolved (main sequence) 
donors we assume that the RLOF will always lead to the 
component merger, thus terminating further evolution 
and aborting potential BH-BH formation. For evolved 
donors, we check whether the RLOF may lead to dynam- 
ical instability. If the donor mass is larger than twice the 
accretor mass {qcrit = 2), we apply standard common en- 
velope prescription l| Webbinkll 1 984^1 with 100% efficiency 
of inspiral orbital energy conversion into the envelope 
ejection. Otherwise, we use the non-conservative evo- 
lution with half of the transferred material leaving the 
system with the angular morn entum specific to a given 
binary ijBelczvnski et al.]l2002() . In the case of accretion, 
the unevolved stars are rejuvenated and may reach larger 
radii than their initial mass would have suggested, while 
the evolved stars are allowed only to increase their mass. 
Once a BH-BH binary is formed, the orbit dec ay time 
due to a gravitational wave emission is calculated ijPetersI 

3. RESULTS 

3.1. Formation and Properties of BH-BH Population 

A large set (iVtot = 10^) of Population HI binaries, 
described in the previous section, is evolved. Evolution 
leads to efficient formation of BH-BH systems (38%), de- 
spite the fact that significant fraction of binaries cease 
to exist when one of the components is disrupted in a 
pair instability SN (46%), or components merge in RLOF 
event (16%). 

Major evolutionary channels along with BH-BH for- 
mation efficiencies are listed in Table 1. Most of the pri- 
mordial systems form on wide orbits and never interact 
(channel bhbhOl). Systems formed on the tight orbits 
interact twice, as first the primary then the secondary 
overfill their Roche lobes (bhbh02, bhbhOS). Systems 
formed on the intermediate orbits interact only once; de- 
pending on the maximum component radii and the mass 
ratio RLOF is initiated cither by the secondary (bhbh04) 
or the primary (within bhbh05). 

In Figure 1 we present distribution of the coalescence 
times for the BH-BH populations. A significant frac- 
tion of BH-BH systems (0.17) coalesces within the Hub- 
ble time (15 Gyr). Binaries with short and intermedi- 
ate initial periods interact and tend to either form tight 
BH-BH or merge in RLOF. Various subpopulations of 
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Fig. 1. — The distribution of BH-BH coalescence times (thick 
solid line) normalized to unity. Additionally, main subpopulations, 
which have evolved to BH-BH stage along different evolutionary 
channels, are presented with thin lines; bhbhOl long dashed dotted, 
bhbh02 dotted, bhbh03 solid, bhbh04 short dashed. Note that 
the subpopulation of systems formed along bhbh05 channels with 
log(tcoal) ~ 3 — 10 is not shown. For definition of channels see 
Table 1. 



BH-BH systems are shown in Fig 1. The stronger (CE 
in case of bhbh02 as compared to non-conservative mass 
transfer) or more frequent interactions (two interactions 
in bhbhOS versus one in bhbh04) lead to shorter coales- 
cence time of BH-BH binary. Large coalescence times 
are found for systems which never interacted, and basi- 
cally remained on the unchanged wide orbits throughout 
the evolution (bhbhOl). Systems with large coalescence 
times dominate the population, since they have formed 
without interactions (wide orbits) and avoided possibility 
of component merger in RLOF events. 

Binary BHs found in our calculation cover a wide 
range of masses; 50-500 and 40-620 Mq for the first 
and the second BH formed in a system, respectively. 
Total mass of BH-BH binaries spans the wide range 
M — 100 — 1000 M0, with most of the systems form- 
ing with M = 100 — 200 M©. The distribution of masses 
of the remaining population peaks at M — S50 M© and 
has a tail extending up to M = 1000 Mq. 

3.2. Observability in Gravitational Waves 

A coalescence of two black holes proceeds through 
three consecutive phases: inspiral, merger, and ring- 
down. The signal to noise ratio {S/N) from each phase 
of a coalescence in the existing a nd near-future inter- 
ferom eters has been calculated by ({Flanagan fc Hughed 
119981) . The typical total mass of a black hole binary in 
the population considered in this Letter is quite large and 
reaches above a few hundred Mq. The inspiral signals 
in interferometric detectors for such massive binaries is 
very low since it falls outside the maximum sensitivity 
range. However, the merger and ringdown signals for 
the advanced LIGO peak for the redshifted total masses 



TABLE 1 0.08 

;h formation channels 



Evolutionary Sequence'' Efficiency 

BHl BH2 0.67 Q 06 

MTl BHl CE2 BH2 0.15 
MTl BHl MT2 BH2 0.07 
BHl MT2 BH2 0.05 8 

all others 0.06 C*i 
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Fig. 2. — The expected observed coalescence rate of Population 
III BH-BH binaries as a function of signal to noise ratio for ad- 
vanced LIGO detector. The thick lines correspond to detection of 
the ringdown signal and the thin lines to detection of the merger 
(burst) signal. Solid lines are for isolated orbital decay of BH-BH 
systems, while dashed lines represent the case when binaries are 
hardened in dense environment (see text for details). 

between 100 and 2000 Mq. 

We calculate the observ ed BH-BH coalesce nce rates 
following the formalism in iBulik "eTall l|200^ . see also 
iBulik & Bclczvnski (20Q3). We estimate the comovmff 
rate of Population HI star formation i?sfr which would 
satisfy the following assumptions: i) formation of Pop- 
ulation HI stars occurs at a constant rate at redshifts 
z from 30 to 10, ii) /mass = 10"'^ of baryonic mat- 
ter g oes into the Population HI stars (Madau & Rccs 
I2OOIV We adopt a flat cosmology model with density 
parameter of matter = 0.3, density parameter of 
cosmological constant JIa = 0.7 and for Hubble constant 
Ho = 100/ikm s-iMpc"\ with h = 0.65. We obtain 
i?sfr — 1-4 X 10~^MQMpc^'^yr^^. Furthermore, we as- 
sume that the binary fraction in Population III stars is 
fb = 0.1, and that the IMF below 100 Mg is flat, and 
extends down to IM©. We flrst calculate the differen- 
tial merger rate as a function of redshift df coai{z) / dAI 
taking into account the delay between formation and co- 
alescence due to gravitational inspiral. The differential 
coalescence rate per unit observed mass is 



dR 

dMo, 



dfcoal{z) 1 dV 

dM 1 + z dz 



dz 



(1) 



where Mobs = M(l -I- z) is the observed (redshifted) to- 
tal mass, zm is the maximum redshift out to which a 
binary is observable, and dV/dz is the comoving volume 
element. The maximum red shift zm is estimated using 
the S/N values estimated bv lFlanagan fc HughesI l|1998f ) 
for the advanced LIGO detector for merger and ringdown 
phases. 

We present the expected rates for advanced LIGO as a 
function of the value of the S/N threshold as solid lines 
in Figure |21 The curves calculated assuming detection 
in merger and ringdown phases are similar because of 



similarity of the dependence of S/N one the masses of 
the system in the two cases. However, it has to be noted 
that the merger S/N is rather uncertain since its calcula- 
tion requires precise prediction of detector noise curve as 
well as knowledge of highly uncertain binary parameters 
(e.g., BH spins). Thus the detection efflciency calculated 
with the use of the predicted merger S/N may be much 
smaller than presented in Figure |21 On the other hand 
the ringdown signal is much better constrained and easier 
to predict. 

Population III BH-BH binaries might have populated 
dense stellar environments in young galaxies. If this was 
the case their coalescence times could have been signifi- 
cantly shortened due to the additional orbital decay by 
three-body interactions. We estimate the observed coa- 
lescence rates in the alternative model in which we allow 
for additional orbital decay due to the dynamical hard- 
ening through simple orbital shrinkage for all BH-BH 
binaries. Orbits are tightened by factor of 10, causing 
the decrease of the coalescence times by factors ~ 10^. 
The corresponding rates are shown in Figure [3 Harden- 
ing does not have a strong impact on the observed rates 
because of two opposing effects. On one hand the tight 
binaries merge at higher redshifts and their detectabil- 
ity drops. On the other, the large population of wide 
non coalescing (for standard model, see Fig. ^ binaries 
is shifted to shorter coalescence times, and in particu- 
lar some may add to the predicted detection rate when 
hardening is included. 

Redshifted total mass of the system is the principal 
quantity that can be inferred from the observation of a 
ringdown signals. In Figure U] we present the observed 
rate as a function of the redshifted total mass, requiring 
a detection with S/N =10. The typical redshifted total 
mass lies between 600 and 1000 Mq. In the case of the 
alternative model with hardening included the typical 
value shifts down to « 300 Mq. This is due to the fact 
that originally long lived binaries contain lighter black 
holes. A calculation using the merger signals leads to 
very similar results. 

4. DISCUSSION 

The principal result of the calculations is presented in 
Figure|21 The expected coalescence rate of Population III 
intermediate-mass BH-BH binaries is high. We predict 
that advanced LIGO should observe above a thousand 
strong {S/N > 10) events per year. What are the uncer- 
tainties of this prediction? 

In our calculation we have assumed that Population 
HI stars were formed at redshifts z from 30 to 10 out 
of /mass = 10"'^ of the baryou mass contained in the 
Universe and that the binary fraction of the initial pop- 
ulation was fb = 0.1. The exact length and duration of 
the Population III star formation era does not affect the 
rate calculation. However, the rate scales linearly with 
/mass and fb, and although we have chosen rather con- 
servative values for these two quantities, the rate may 
decrease if the adopted values were significantly lowered. 

The IMF of Population III stars is another unknown. 
Numerical investigations show that the IMF leans to- 
wards massive stars. We assumed rather steep IMF (with 
the slope a — —2), to assure that we do not overproduce 
highest mass, and therefore the easiest to detect, BH- 
BH binaries. As it turns out, the change of IMF slope 
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Fig. 3. — The difTcrcntial rate as a function of observed (red- 
shifted) total mass of BH-BH system. In this calculation we re- 
quired a detection with high signal to noise ratio {S/N > 10). 

does not significantly alter the detection rate. Farther 
steepening of the IMF (a = —3) decreases the expected 
rate by half, while flatter IMF assumption (a = —1) in- 
creases the expected rate by a factor of two. Note that 
we maintain a flat IMF between 1 M© and 100 Mq when 
calculating fsim ~ the fraction of stars that we simulate 
out of the total population. The observed rate scales lin- 
early with fsim- We also fixed the maximum mass of the 
Population III stars to 500 M©. Had we allowed for pos- 
sibility of star formation with the higher mass the rate 
would increase. 

The estimate of the lifetime of the black hole bina- 
ries can be strongly affected if the binaries are hardened 
by interactions in dense stellar environments. Our sim- 
ulations show that a large number of black hole bina- 
ries should be formed with the lifetimes in excess of the 
Hubble time. Interactions in dense systems may signif- 
icantly shorten their lifetimes. We demonstrated above 
that hardening does not affect strongly the expected rate. 
Another possibility arises that the interactions disrupt 
some of BH-BH binaries. This would deplete the popu- 
lation of wide systems. None of the two effects, unless 
operating on extremely short timescales, can affect much 
systems with shortest coalescence timescales. Therefore, 
a combination of the two effects may constrain the merg- 
ers of Population III BH-BH binaries to large redshifts. 



Abt, H. A. 1983, ARA&A, 21, 343 

Baraffe, I., Heger, A., & Woosley, S. E. 2001, ApJ, 550, 890 
Belczynski, K., Kalogera, V., & Bulik, T. 2002, ApJ, 572, 407 
Bromm, V., Coppi, P. S., & Larson, R. B. 1999, ApJ, 527, L5 
— . 2002, ApJ, 564, 23 

BuUk, T. & Belczynski, K. 2003, ApJ, 589, L37 

BuUk, T., Belczynski, K., & Rudak, B. 2004, A&A, 415, 407 

Chabrier, G. 2003, PASP, 115, 763 

Flanagan, E. E. & Hughes, S. A. 1998, Phys. Rev. D, 57, 4535 
Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, 

D. H. 2003, ApJ, 591, 288 
Heger, A. & Woosley, S. E. 2002, ApJ, 567, 532 



In order to estimate the S/N for advanced L IGO we 
have used the formulae of iFlanagan fc HughesI (^QS). 
The predicted values of S/N may still change when more 
realistic noise curves and binary gravitational wave form 
templates are known. Because of large masses of the 
binaries considered in this work the changes in the low- 
frequency range are most important. The typical ring- 
down frequency is Vqnr ~ 9O(3OOM0/M) Hz, so the rate 
is most sensitive to the detector performance in the low 
frequency region. The scaling of the rate with the change 
of S/N normalization can be read off Figure |2| 

The influence of the evolutionary model assumptions 
on detection rate was tested. We changed the CE effi- 
ciency (increase/decrease by factor of 2); altered evolu- 
tion through stable RLOF phases from fully conserva- 
tive to non-conservative cases; changed the specific an- 
gular momentum of the matter leaving the system dur- 
ing RLOF (increase/decrease by factor of 2) and var- 
ied the critical mass ratio over which the dynamical in- 
stability develops (from qcrit = 2 in standard model to 
9crit = 1 — 3). The detection rate was decreased at most 
by a factor of 3 in the above models. Therefore, the 
predicted rate does not depend strongly on the binary 
evolution within the model assumptions. 

Summarizing, in our calculations we have tried to use 
rather conservative assumptions in order not to overesti- 
mate the detection rate of Population HI BH-BH merg- 
ers. If several of the assumptions and values of the model 
parameters are changed within reasonable limits, we still 
obtain a significant detection rate. Even if the rate pre- 
dicted for our already conservative model is decreased by 
2-3 orders of magnitude to allow for different aforemen- 
tioned uncertainties, we still are left with several strong 
events per year for advanced LIGO detector. For initial 
LIGO phase the detection rate falls below one event per 
year. 

We have shown that Population III stars may lead 
to formation of a large number of binaries containing 
intermediate-mass black holes. A significant fraction of 
such systems has coalescence times smaller than the Hub- 
ble time. Coalescences of intermediate-mass BH-BH bi- 
naries should be detectable by advanced interferomctric 
gravitational wave detectors during the ringdown phase, 
and possibly also merger phase, provided that accurate 
templates are available. Given the large expected rate of 
observed coalescences such events could be the primary 
targets for LIGO burst search. 
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